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Abstract 

Results of torsional oscillation tests are reported that were performed at the tempera- 
ture T = 230 °C on melts of a hybrid nanocomposite consisting of isotactic polypropylene 
reinforced with 5 wt.% of montmorillonite clay. Prior to mechanical testing, specimens 
were annealed at temperatures ranging from T a = 250 to 310 °C for various amounts 
of time (from 15 to 420 min). Thermal treatment induced degradation of the matrix 
and a pronounced decrease in its molecular weight. An integro-differential equation is 
derived for the evolution of molecular weight based on the fragmentation-aggregation 
concept. This relation involves two adjustable parameters that are found by fitting ob- 
servations. With reference to the theory of transient networks, constitutive equations are 
\ developed for the viscoelastic response of nanocomposite melts. The stress-strain rela- 

tions are characterized by three material constants (the shear modulus, the average energy 
for rearrangement of strands and the standard deviation of activation energies) that are 
determined by matching the dependencies of storage and loss moduli on frequency of oscil- 
lations. Good agreement is demonstrated between the experimental data and the results 
of numerical simulation. It is revealed that the average energy for separation of strands 
from temporary junctions is independent of molecular weight, whereas the elastic modu- 
lus and the standard deviation of activation energies linearly increase with mass-average 
molecular weight 
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Introduction 



This paper is concerned with (i) the kinetics of thermal degradation of a hybrid nanocomposite 
with isotactic polypropylene (iPP) matrix reinforced with montmorillonite (MMT) clay and (ii) 
the effect of degradation on the viscoelastic response of the nanocomposite in the melt state. 

The choice of iPP for the experimental analysis is explained by numerous industrial applica- 
tions of this semicrystalline polymer (oriented films for packaging, reinforcing fibres, non-woven 
fabrics, pipes, etc.). Montmorillonite is an inorganic clay conventionally used for preparation of 
hybrid nanocomposites. It possesses a layered structure constructed of two tetrahedral sheets 
of silica surrounding an octahedral sheet of alumina or magnesia. The layers (with a thickness 
of 1 nm) are stacked by weak dipole forces, while the galleries between the layers are occupied 
by metal cations. 

The focus on the influence of thermal degradation of a polymer-clay nanocomposite on its 
time-dependent response may be explained by two reasons: (i) from the standpoint of appli- 
cations, the effect of degradation on the viscoelastic behavior is of essential importance for 
the prediction of mechanical properties of reprocessed industrial and post-consumer plastic 
wastes pQ; (ii) from the point of view of fundamental research, thermal degradation of a poly- 
mer is tantamount to a reduction in its molecular weight [2], which implies that the analysis 
of the time-dependent response of nanocomposite melts annealed at various temperatures T a 
for various amounts of time t a sheds some light on their structure-property relations (correla- 
tions between the molecular weight of the matrix and the material parameters describing the 
mechanical response). 

The kinetics of thermal and thermo-oxidative degradation of polypropylene has attracted 
substantial attention in the past decade, see [3 EJ E] , to mention a few. Annealing, thermal 
degradation and stability of nanocomposites with a polypropylene matrix and MMT filler were 
recently studied in [7J |Hl HOI HH U21 UHj • Some relations between the viscoelastic response of 
iPP-MMT nanocomposites and their structure have been established in |14 [ IT5 ^ I17j. 

Despite substantial progress in our understanding of the degradation process and its effect 
on the mechanical response of hybrid nanocomposites, it is difficult, however, to mention a con- 
stitutive model that adequately predicts the time- dependent behavior of iPP-MMT nanocom- 
posite annealed for a given time t a at a required temperature T a above the melting temperature 
T 

The objective of this study is three-fold: 

1. To report experimental data in isothermal torsional oscillation tests with small strains 
(at the temperature T = 230 °C, which is a typical temperature for injection-molding of 
iPP) on specimens annealed for various amounts of time t a (ranging from 15 to 420 min) 
at various temperatures T a (in the interval from 250 to 310 °C). 

2. To develop kinetic equations for the evolution of number-average and mass-average molec- 
ular weights at thermal degradation and to find adjustable parameters in these relations 
by fitting the experimental data. 

3. To derive constitutive equations for the viscoelastic behavior of a nanocomposite melt, 
to determine material constants in the stress-strain relations by matching the dependen- 
cies of the storage and loss moduli on frequency, and to evaluate the effect of thermal 
degradation on the time-dependent response. 
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To make the model tractable from the mathematical standpoint, we adopt the homogeniza- 
tion hypothesis. According to it, a complicated micro-structure of a nanocomposite may be 
replaced by an equivalent phase, whose response captures essential features of the mechanical 
behavior of the nanocomposite. As the equivalent phase, a network of macromolecules is chosen 
in the present study. 

Following ^Hl 023 120] , degradation of the host matrix in an hybrid nanocomposite is treated 
as a combination of two thermally-activated processes: (i) binary scission of chains, and (ii) 
annihilation of end- and side-groups (diffusion of small-size fragmentation products and their 
subsequent evaporation through the surface of a specimen). To simplify the analysis, we accept 
the conventional assumptions that (i) the probability of a scission event is independent of a 
chain's length and the position of a bond along the backbone of a chain; (ii) the diffusivity of 
separated end- and side-groups is so large at the exposure temperature T a that the kinetics of 
diffusion of oligomers and their evaporation may be disregarded. As a result, the number of 
material constants in the kinetic equation is reduced to two. These quantities are found by 
matching a master-curve for the decrease in the mass-average molecular weight with exposure 
time t a . 

To develop constitutive equations for a nanocomposite melt, we accept the concept of tran- 
sient networks fZT\ I22J G31 121] • The melt is treated as a network of strands bridged by temporary 
junctions (entanglements and physical cross-links on the surfaces of clay platelets). At random 
times, active strands separate from their junctions being excited by thermal fluctuations, and 
dangling strands merge with the network. Following [23 |2E|, we assume the network to be 
strongly heterogeneous in the sense that different junctions are characterized by different acti- 
vation energies for detachment of strands. The theory of temporary networks has been applied 
to describe the time-dependent behavior of polypropylene in the solid state in (23 I2H EHl EE] , 
to mention a few. The previous studies, however, did not pay much attention to the effect 
of molecular weight on material parameters in the constitutive equations. The fact that the 
viscoelastic response of polypropylene melts is strongly affected by the distribution of chains' 
lengths has been established experimentally in [23 HUB EH] E21 ESU Ell- The aim of this work 
is to evaluate the influence of mass-average molecular weight on (i) the concentration of active 
strands in an equivalent network, and (ii) the distribution function for temporary junctions 
with various activation energies. 

The exposition is organized as follows. First, experimental data in torsional oscillation tests 
are reported on annealed specimens, and their mass-average molecular weight is determined 
by using observations for complex viscosity. Afterwards, kinetic equations are developed for 
changes in the concentration of chains with various lengths induced by thermal treatment. 
Adjustable parameters in these equations are determined by fitting the observations for molec- 
ular weight. We proceed with the derivation of stress-strain relations for an heterogeneous 
transient network at three-dimensional deformations with small strains. The constitutive equa- 
tions involve three material constants that are found by matching the data for storage and loss 
moduli as functions of frequency of oscillations. Finally, we establish correlations between the 
adjustable parameters in the stress-strain relations and the mass-average molecular weight of 
the polypropylene matrix and discuss the physical meaning of these relationships. 
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Experimental procedure 



Isotactic polypropylene PP 1012 (density 0.906 g/cm 3 , melt flow rate 1.2 g/10 min) was pur- 
chased from BP Amoco Polymers, Inc. Concentrate C-44PA containing 43 wt.% of intercalated 
montmorillonite clay in de-agglomerated form was supplied by Nanocor Co. 

To prepare hybrid nanocomposite with 5 wt.% of MMT clay, appropriate amounts of iso- 
tactic polypropylene and the concentrate were dried overnight at the temperature T = 100 °C, 
pellets were mixed and melt-blended in a twin-screw extruder (Brabender Instruments, Inc.) 
with a screw rate of 30 rpm and temperatures in the extruder barrel of 260, 290, 300 and 290 °C 
from hopper to die, respectively. Strands from the extruder were cooled in a water bath, cut 
with a pelletizer, and dried in an oven at 100 °C for 12 h. Circular plates with diameter 64 
mm and thickness 3 mm were molded in injection-molding machine Battenfeld 1000/315 CDC 
(Battenfeld). Specimens for mechanical tests (with diameter 30 mm) were cut from the plates. 

Our choice of the concentration of MMT clay in the hybrid nanocomposite (0 = 5 wt.%) 
is explained by the fact that this amount of nanoclay is sufficient to improve substantially 
mechanical properties of neat polypropylene. According to [IB], the tensile strength of iPP 
MMT nanocomposite reaches its maximum when the clay concentration equals 5 wt.%. Another 
reason for this choice is that reinforcement of conventional polymers with MMT clay at higher 
concentrations of filler practically does not improve their thermal stability |3*3] . 

To evaluate the melting temperature T m of the matrix and the nanocomposite, DSC (differ- 
ential scanning calorimetry) measurements were preformed by using DSC 910S apparatus (TA 
Instruments). The calorimeter was calibrated with indium as a standard. Two specimens of 
neat iPP and iPP-MMT nanocomposite with weights of approximately 13 mg were tested with 
a heating rate of 10 K/min from room temperature to 200 °C under nitrogen. The melting 
temperature T m = 172 °C was determined for isotactic polypropylene as the point correspond- 
ing to the peak on the melting curves. No substantial changes in the melting temperature 
were found for the nanocomposite. This conclusion is in agreement with the results of previous 
studies [13J. 

To analyze lattice spacing in the montmorillonite clay and changes in the crystalline mor- 
phology of iPP driven by the presence of filler, X-ray diffraction tests were performed on iso- 
tactic polypropylene and polypropylene-clay nanocomposite by using Rugaku D-max B diffrac- 
tometer with Cu-K a radiation (A = 1.54 A) generated by a tube with a voltage of 40 kV and 
a current of 30 mA. The Bragg scattering angle ranged from 20 = 3 to 20 = 60° with the step 
of 0.06°. The diffraction spectrum of the hybrid nanocomposite revealed two peaks: a rather 
large and wide intensity maximum at 20 = 3.0° (the basal spacing of 30 A) and a smaller 
one at 20 = 6.5° (the basal spacing of 17 A). The presence of these peaks indicated a high 
level of intercalation of polymer chains between clay platelets. These observations are similar 
to those found by other researchers |T4^ ITo! IT7| I3*o] on iPP-MMT nanocomposites with similar 
concentrations of clay. 

Among other interesting features of the X-ray diffractograms, we would mention: (i) an 
increase in a(110) peak of the nanocomposite compared to that of iPP (by a factor of 5), (ii) 
the growth of a(040) peak (by a factor of 3.3), (iii) a decrease in a(130) peak (by 26 %), (iv) 
the formation of 7(117) peak at 20 = 19.6°, (v) the disappearance of a(lll) and a{— 131) 
peaks, and (vi) a strong growth of a(060) peak at 20 = 25.2° (by a factor of 3.5). 

It is worth noting that a pronounced reduction in the diffraction peaks at 20 = 21.2° 
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was recently reported in jSTj for hot-stretched polypropylene samples. It was associated with 
preferential orientation of crystallites driven by deformation of specimens. With reference to this 
assertion, the fact that the a(lll) and a(— 131) peaks were not observed in the nanocomposite 
may be ascribed to the effect of anisotropically distributed clay platelets on the orientation of 
lamellar blocks in injection-molded samples. 

Dynamic tests were performed by using RMS-800 rheometric mechanical spectrometer with 
parallel disks (diameter 25 mm, gap length 2 mm). The shear storage modulus G' and the 
shear loss modulus G" were measured in oscillation tests (the frequency-sweep mode) with the 
amplitude of 15 % and various frequencies uj ranging from 0.1 to 100 rad/s. The choice of the 
amplitude of oscillations was driven by the following requirements: (i) mechanical tests were 
performed in the region of linear viscoelasticity, and (ii) the torque was less than its ultimate 
value 0.2 N-m. The limitation on the minimum frequency of oscillations was imposed by the 
condition that the torque exceeded its minimum value 2.0 • 10 -4 N-m. To check that the storage 
and loss moduli were not affected by the strain amplitude, several tests were repeated with the 
amplitude of 5 %. No change in the dynamic moduli was observed. The temperature in the 
chamber was controlled with a standard thermocouple that showed that the temperature of 
specimens remained practically constant (with the accuracy of ±0.5 °C). 

Prior to mechanical tests, specimens were annealed in the spectrometer at the temperatures 
T a = 250, 270, 290 and 310 °C (with gap length 3 mm) for various amounts of time t a ranging 
from 15 min to 7 h (420 min). After thermal treatment of a specimen, the temperature was 
reduced to the test temperature T = 230 °C, the specimen was thermally equilibrated (during 
5 min), the gap length was reduced to 2 mm, an extraneous material was carefully removed, 
and the shear storage and loss moduli were measured at various frequencies uj starting from 
the lowest one. Each test was performed on a new specimen. 

We suppose that squeezing of samples between plates of the spectrometer and removal 
of the extraneous material substantially reduced the effect of thermo-oxidative degradation 
(compared to that of thermal degradation) on the mechanical response, because the major part 
of the material where oxidative degradation occurred was taken away (about one third of the 
initial mass of each specimen). However, we cannot exclude entirely the effect of diffusion of 
oxygen to the central part of the samples, in particular, at the highest temperatures (T a = 290 
and 310 °C) used in the experiments [38J. 

The storage G' and loss G" moduli, as well as the modulus of the complex viscosity r\ are 
plotted versus the logarithm (log = log 10 ) of frequency uj in Figures 1 to 12 (conventional semi- 
logarithmic plots are used to characterize changes in these quantities with frequency). The 
shapes of the curves presented coincide qualitatively with those reported by other researchers, 
see, e.g., [T5] . Given an annealing time t a and an annealing temperature T a , the storage modulus 
G' and the loss modulus G" strongly increase with frequency, whereas the complex viscosity r\ 
noticeably decreases with uj. For a fixed frequency uj, the dynamic moduli are pronouncedly 
reduced with t a and T a . 

Evaluation of molecular weight 

To assess changes in the molecular weight of the hybrid nanocomposite induced by thermal 
degradation of the matrix, the experimental dependencies T](uj) depicted in Figures 3, 6, 9 and 
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12 are approximated by the Cross model 



Ar] 

id*) = voo + YTJ^r^ (1) 

where 77^ is the high-frequency complex viscosity (u — > 00), Ar] = r] — r]^, r] is the zero- 
frequency complex viscosity (u — 0), and a and r are adjustable parameters. 

Each curve r](u) is approximated separately. To find the constants rjoo, Ar], a and r in Eq. 
(1), we fix some intervals [0, a max ] and [0, r max ], where the "best-fit" parameters a and r are 
assumed to be located, and divide these intervals into J subintervals by the points «W = iAa 
and t^' = j At (i, j = 1, . . . , J — 1) with Aa = a max /J and Ar = r max /J. For any pair 
{a''',T^}, the coefficients and Ar/ in Eq. (1) are found by the least-squares method from 
the condition of minimum of the function 



where the sum is calculated over all experimental points oj m depicted in Figures 3, 6, 9 and 
12, ?7 ex p is the complex viscosity measured in a test, and 77 num is given by Eq. (1). The "best- 
fit" parameters a and r are determined from the condition of minimum of the function F on 
the set {a^\ = 1, . . . , J — 1)}. After finding the "best-fit" values a' 8 ' and this 

procedure is repeated twice for the new intervals and , to ensure an 

acceptable accuracy of fitting. Figures 3, 6, 9 and 12 demonstrate good agreement between the 
experimental data and the results of numerical simulation. 

After finding the zero-frequency viscosity 770, the mass- average molecular weight M w is 
determined by the conventional equation 

3.4 

(2) 



rjo_ = ( M w \ 



where 77g ef and M^ ef are the zero-frequency viscosity and the mass-average molecular weight of a 
reference (not subjected to thermal treatment) specimen. The ratio of mass-average molecular 
weights 

Mw 
w M ref 

w 

is found from Eq. (2) by using the experimental data for r/ - This method of determining the 
ratio of mass-average molecular weights of polypropylene was previously used in 40 . 

It is worth noting that Eq. (2) is traditionally employed for the evaluation of molecular 
weight of neat polymers. Its applicability to melts of hybrid nanocomposites is grounded on 
the conventional models in rheology of particulate suspensions [H], which presume that the 
relative viscosity (the ratio of the viscosity of a suspension to that of the neat polymer melt) is 
independent of the melt's structure and is determined by the concentration of filler exclusively. 

The ratio d w is plotted versus annealing time t a in Figure 13. With reference to jHl l2"Uj . 
we accept the time-temperature superposition principle for degradation of the nanocomposite. 
According to this hypothesis, the dependencies d w (t a ) measured at various annealing temper- 
atures T a and plotted in semi- logarithmic coordinates (d w versus logt a ) may be superposed 
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with an acceptable level of accuracy by shifting the observations along the time-axis. Apply- 
ing this approach, we construct the master-curve depicted in Figure 13. The experimental 
data at T^ ef = 290 °C are presented without changes. Observations at the other temperatures 
(T a = 250, 270 and 310 °C) are shifted along the time-axis by appropriate amounts A that are 
determined from the condition that the experimental data produce a smooth master-curve. 

The parameter A is plotted versus the absolute temperature T a in Figure 14. The experi- 
mental data are approximated by the Arrhenius dependence 

\nA = A -^, (3) 

where the coefficients A i (i = 0, 1) are determined by the least-squares method. Figure 14 
demonstrates that Eq. (3) ensures quite acceptable fit of the observations. 

For each annealing time t a and annealing temperature T a , we (i) calculate the material 
constants a and r in Eq. (1) by matching the observations depicted in Figures 3, 6, 9 and 12, 
(ii) find the ratio of mass-average molecular weights d w from Eqs. (1) and (2), and (iii) plot 
the adjustable parameters a and r versus d w in Figures 15 and 16. The experimental data are 
approximated by the phenomenological relations 

a = a — axd w , r = t + r ± d w , (4) 

where the coefficients ctj and r, (i = 0, 1) are found by the least-squares technique. Figures 
15 and 16 reveal that Eq. (4) provides reasonable quality of matching the observations. The 
exponent a in Eq. (1) slightly decreases with mass- average molecular weight, whereas the 
characteristic time r strongly grows with M w . 



Kinetic equations 

Our aim now is to develop kinetic equations for thermal degradation of a hybrid nanocomposite 
and to find adjustable parameters in these relations by matching the experimental data depicted 
in Figure 13. For this purpose, we (i) adopt a homogenization method, according to which a 
complicated micro-structure of a polymer-clay nanocomposite may be replaced by an equiva- 
lent network of macromolecules, and (ii) accept the fragmentation-annihilation concept. The 
latter means that the degradation process is treated as a combination of two thermally-induced 
processes: binary fragmentation of chains and annihilation (subsequent breakage, diffusion and 
evaporation through the surface of a sample) of end- and side-groups. 

Binary scission of chains 

Denote by N(t) the number of macromolecules per unit mass of an equivalent network at an 
arbitrary instant t > 0. Following common practice, we treat chains as sequences of segments 
connected by bonds. Denote by N^(t) is the number of chains (per unit mass) at time t 
containing k segments (k — 1, 2, . . .). The functions N k (t) obey the conservation law 

oo 

#(f) = (5) 
k=i 
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Binary scission (fragmentation) of chains is described by the reactions 

Nk^Ni + N^ (Z = 1, — 1). 

Denote by 7 the rate of scission (the number of scission events per bond between segments 
per unit time). Assuming 7 to be a function of temperature T only (which implies that 7 is 
independent of the number of segments in a chain), we arrive at the kinetic equations for the 
functions N k (t) 

HIV 00 

a -±(t) = - 7 (* - l)N k (t) + 2 7 £ iV,(t). (6) 

j=k+l 

The coefficient k — 1 in the first term describes the number of possible scission events in a chain 
containing k segments. The coefficient "2" before the sum in Eq. (6) indicates that there are 
two opportunities ("left" and "right") to obtain a chain with k segments after scission of a 
chain with a larger number of segments. 



Annihilation of chains 

Thermal fluctuations in a network induce not only binary scission of macromolecules, but also 
detachment of end- and side-groups from polymer chains. As these groups are rather small, 
they have relatively large diffusivity, and can easily leave a polymer specimen. A decrease 
in a sample's mass with time driven by separation and subsequent desorption of end- and 
side-groups and is treated as their annihilation. 

Following [HHH], we suppose that detachment of small groups within the interval [t,t + dt] 
may be thought of as transformation of a chain with k segments into a chain with k — 1 
segment. Denote by Tk the ratio of the number of chains with k segments that lose a segment 
per unit time (due to the annihilation process) to the entire number of these chains N k . The 
parameter T k is proportional to the number of thermal fluctuations (per unit time) that induce 
detachment of side-groups and the number of side-groups per macromolecule. As both these 
quantities linearly increase with k, one can write 

T k = Tk 2 , (7) 

where T is a temperature-dependent material parameter. The kinetic equation for the fragment- 
ation-annihilation process reads 

dN 00 

< ^±{t) = - 1 {k-l)N k {t)+2 1 J2 N,(t) + T[(k + l) 2 N k+1 (t)-k 2 N k (t) 

j=k+l 

Introducing the concentrations of chains with k segments, 

N k {t) 



n k (t) 



N 



where Nq = N(0) is the total number of chains at the initial instant t = 0, we present this 
equation in the form 



dn k 

^ j=k+i 



(t) = - 7 (fc - l)n k (t) + 2 7 £ n^ + TUk + lfn^^-en^t)]. (9) 
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An explicit solution 

Our purpose now is to analyze changes in the number- average molecular weight M n and the 
mass-average molecular weight M w determined by the conventional relations 

M (t) = ^°=i krik ^ M (t)(t) = ^fc°=i fMQ ( 10) 

when the functions nk(t) are governed by Eq. (9) with an arbitrary initial condition 

n k (0)=n ok (k = 1,2, . . .). 

Following common practice, it is convenient to suppose that the number of segments in a chain 
is large compared to unity and to replace the discrete index k in Eq. (9) by a continuous 
argument x. This results in the integro-differential equation for the function n(t,x), 

On f°° / \ 

— (t,x) = -'yxn(t,x)+2'yj n(t,y)dy + T—[x 2 n(t,x)), ra(0, x) = n (x), (11) 

where Uq(x) is a given function. We do not formulate boundary conditions for the function 
n(t,x), but assume that this function does not grow very strongly at x = and decays rapidly 
at x — > oo in the sense that the integrals exist 

M m (t)= x m n{t,x)dx (m = 0,1,2,...). (12) 
Jo 

Multiplying Eq. (11) by x m (m = 0, 1, . . .), integrating over x, and using notation (12), we find 
that 



it :(t) = - 1 M m+1 (t) + 2 1 J o x m dxj x n(t,y)dy + Tj o x rn —(x 2 n{t,x))dx. (13) 
The first integral is transformed by changing the order of integration, 

roc roo poo py | poo | 

/ x m dx / n{t,y)dy= n{t,y)dy x m dx = — / y m+1 n{t,y)dy = —M m+1 (t). 

Jo Jx Jo Jo m + 1 Jo m+1 

The other integral is calculated by integration by parts, 

roo d / \ r°° 

/ x m — (x 2 n(t,x))dx = -m / x m+1 n(t,x)dx = -mM m+ At). 
Jo ox ^ 1 Jo 

Substitution of these expressions into Eq. (13) results in the differential equation 

dM m . . (m — 1 \ „ , , N 

it) = -7 — rr + «™ M m+i(*) ( 14 ) 
\m + 1 / 



with 

« = -. (15) 

7 

Assuming the equivalent network of chains to be monodisperse at the initial instant t — 0, 

n (x) = 5(x — L), (16) 
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where L is the initial length of chains, we solve Eq. (14) by Charlesby's method [12] • As 
Eq. (14) is linear with respect to the unknown function n(t,x), appropriate formulas for the 
moments M m (t) corresponding to an arbitrary initial condition n (x) are developed by the 
superposition method. 

The mth moment M m (t) is expanded into the Taylor series in time, 



fc=0 



k\ 



(17) 



where M^(0) stands for the fcth derivative at the point t = 0. Substitution of Eq. (17) into 
Eq. (14) implies that 



fc-i 



M«(0) = (-7) fe II 

j=0 



m + j — 1 
m + j + 1 



+ «(m + j) 



M m+k (0). 



It follows from Eqs. (12) and (16) that 

M m (0) = L m . 

Substitution of these expressions into Eq. (17) results in 



M m (t) = U 



1 + £ Ank(-7tLY 
k=l 



where 



fc-1 r 



A ™k - TJ n 

K - 3=0 



m + j 



+ «;(m + j) 



.171 + j + 1 

Introducing the new variable j' — j + 1 and omitting the prime, we obtain 



A m k — I I 

fJiJ 



" Win ■ j - 2 



m + j 



+ n(m + j - 1) 



(19) 



Equation (18) implies that for an arbitrary initial condition no(x), the moments M m (t) are 
given by 



M m (t) 



no{x)x 



i + 5> 



mfc I 



-7tx) 



fe=l 



dx = M m (0) + £ A mfc M m+fc (0)(- 7 t) fc . (20) 



fc=i 



Although Eqs. (19) and (20) provide explicit expressions for the moments M m (t), they are not 
convenient for the numerical analysis, because the series converges slowly. These formulas are 
helpful, however, for the evaluation of changes in M m (t) at small times, jt <S 1. Neglecting 
terms beyond the first order of smallness in Eq. (20) and using Eqs. (15) and (19), we find 
that 

M (t) = M (0) + Afi(0)7t, Mi(t) = Mi(0) - M 2 (0)Tt, 



M 2 (t) = M 2 (0)-M 3 (0)Q + 2T)t. 
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According to these equations, changes in the moments M (t) and M-y{t) are governed by two 
different processes: an increase in M (t) is driven by fragmentation of chains, whereas a decrease 
in Mi(t) is induced by annihilation of end- and side-groups. Introducing the notation 



dn(t) 



M n (t) 



M n (0)' 

and using Eqs. (10) and (12), we obtain 

d n (t) = 1 - (M n (0) 7 + M w (0)r)t, d w (t) 
where M z (t) = M 3 (t)/M 2 (t). 



d w (t) 



M w (t) 
M w (0) 



M w (o)r + M z (o)g + 2r) 



(21) 



t, (22) 



Fitting of observations 

As the series in Eq. (18) converges slowly, we analyze the evolution of the mass- average molec- 
ular weight M w with elapsed time t numerically. For this purpose, we integrate Eq. (9) with 
the initial condition 

nok = hK, (23) 

where Sij denotes the Kronecker delta. Equation (23) corresponds to a monodisperse distri- 
bution of chains in an equivalent network that contain K segments at the initial instant. We 
chose this assumption, because the precise initial distribution of chains in a nanocomposite is 
unknown. Bearing in mind that if the maximal number of segments in a chain equals K at 
t = 0, no chains with higher number of segments can appear at t > due to the fragmentation- 
annihilation process and using Eq. (15), we re- write Eq. (9) as follows: 

d K 

J£(t) = _ 7 (fc_l) nfc (t) + 2 7 J2 n 3 (t) + K[(k + l) 2 n k+ i(t)-k 2 n k (t)] (k = 1, 2, . . . , K). 
at j=k+i 

(24) 

We fix the value K = 100 and integrate Eq. (24) with the step At = 0.1. This, relatively large, 
step is chosen because the fragmentation rate under consideration 7 is quite small (of the order 
oflO' 8 ). 

To find the adjustable parameters 7 and k, we fix some intervals [0, 7 max ] and [0, « max ], where 
the "best-fit" parameters 7 and k are assumed to be located, and divide these intervals into J 
subintervals by the points 7W = iAj, and = jAn — 1, . . . , J — 1) with A7 = 7 max / J 
and An = « max / J. For any pair {7^, k^}, Eq. (24) with initial condition (23) is integrated by 
the Runge-Kutta method. The best-fit parameters 7 and k are determined from the condition 
of minimum of the function 



^ = E[ < C P (*m)-<C m (*m) 



where the sum is calculated over all times t m at which observations are presented in Figure 13, 
d^ xp is the ratio of mass-average molecular weights measured in the tests, and d™ m is given 
by Eqs. (10) and (21). Figure 13 demonstrates fair agreement between the observations on 
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specimens annealed at various temperatures T a and the results of numerical simulation with 
7 = 7.9 • 10" 8 and k = 436.0. 

To ensure the accuracy of numerical simulation, we use three tests. First, we verify that at 
k = 0, the first moment M% remains independent of time [this conclusion follows from Eq. (14) 
with k = 0]. Secondly, we increase K by twice, decrease the rate of fragmentation 7 by twice 
and check that the moments M n (t) and M w (t) remain unchanged. The latter implies that the 
results of numerical analysis are independent of our choice of K = 100. Finally, we perform 
simulation at relatively small times and confirm that the numerical results for the moments 
M m (t) (m = 0, 1,2) coincide with analytical solution (22). 

The rate of fragmentation 7 found by matching observations on iPP-MMT nanocomposite 
is of the same order of magnitude as that determined in [20] for degradation of polystyrene 
at T a = 275 °C. The parameter k is of the order of 10 2 , which implies that the influence of 
annihilation of side-groups on the degradation process is substantial [20J. 

To compare our results of numerical analysis with observations reported by other researchers, 
we recall that Eq. (3) is based on two hypotheses: (i) the rate of fragmentation 7 follows the 
Arrhenius dependence on the annealing temperature T a , 

7 = 7 exp(-— — ), (25) 



RT a 

where E is the activation energy and R is the universal gas constant, and (ii) the ratio k of the 
rates of annihilation and fragmentation is independent of annealing temperature T a . It follows 
from Eq. (25) that the shift factor A 

A ~ ref ' 

7o 

E „ E 

A " = w ~ b; < 26 > 

Calculating the activation energy E from Eq. (26) and Figure 14, we find that E = 61.8 
kJ/mol. This value is rather close to the activation energies for thermal degradation of isotactic 
polypropylene E = 50 to 120 kJ/mol determined in [Hj based on results of thermo-gravimetrical 
tests at low conversion factors. 



is given by Eq. (3) with 



Constitutive equations 

Our aim now is to fit the experimental data for storage and loss moduli of the hybrid nanocom- 
posite annealed at various temperatures T a . For this purpose, we derive constitutive equations 
for the viscoelastic response of a nanocomposite melt at three-dimensional deformations with 
small strains, simplify these equations for steady shear oscillations, find adjustable parameters 
in the stress-strain relations by matching the observations depicted in Figures 1, 2, 4, 5, 7, 8, 10 
and 11, and analyze the effect of mass-average molecular weight on the material constants. Our 
analysis is based on the assumption that the characteristic time for thermal degradation (of 
the order of a few hours) substantially exceeds the characteristic time for relaxation of stresses 
in a nanocomposite melt (of the order of a few seconds), which implies that scission of macro- 
molecules and annihilation of side-groups may be disregarded in the analysis of mechanical 
tests. 
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With reference to the concept of transient networks, a nanocomposite melt is treated as 
an equivalent network of strands bridged by temporary junctions (entanglements and physical 
cross-links whose life-time does not exceed the characteristic time of a mechanical test). A 
strand whose ends are linked to contiguous junctions is treated as an active one. When an end 
of an active strand separates from a junction, the strand is transformed into the dangling state. 
When a free end of a dangling strand captures a nearby junction, the strand returns into the 
active state. Separation of active strands from their junctions and merging of dangling strands 
with the network occur at random times when the strands are excited by thermal fluctuations. 
According to the theory of thermally-activated processes [13] , the rate of detachment of strands 
from temporary junctions $ is governed by the equation 

$ = $o6XP (-^)' (2?) 

where <3>o is the attempt rate (the number of separation events per strand per unit time), &b 
is Boltzmann's constant, T is the absolute temperature, and v > is the activation energy for 
separation of an active strand. The coefficient $o i n Eq. (27) is independent of the activation 
energy v and is determined by the current temperature T only. Confining ourselves to isothermal 
processes at a reference temperature T and introducing the dimensionless activation energy 
v = v/(k^T rei ), we find from Eq. (27) that 

= <5> exp{-v). (28) 

To describe the time-dependent response of a nanocomposite melt, we follow the approach 
proposed in I2H] and suppose that different junctions are characterized by different di- 
mensionless activation energies v. The distribution of active strands in a transient network is 
determined by the number of active strands per unit mass N a and the distribution function p{v) . 
The quantity N a p(v)dv equals the number of active strands per unit mass linked to junctions 
with the dimensionless activation energies u belonging to the interval [v,v + dv]. 

Separation of active strands from temporary junctions and merging of dangling strands with 
the network are entirely described by the function u(t, r, v) that equals the number (per unit 
mass) of active strands at time t > linked to temporary junctions with activation energy v 
which have last merged with the network before instant r e [0,t]. 

The quantity u(t, t, v) equals the number of active strands (per unit mass) with the activation 
energy v at time t, 

v(t,t,v) = N a p(v). (29) 



The function 

. . dv . 



(30) 



determines the rate of reformation for dangling chains: the amount (p(r, v)dr equals the number 
of dangling strands (per unit mass) that merge with temporary junctions with activation energy 
v within the interval [r, r + dr]. The quantity 

dv 

—(t,r,v) dr 
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is the number of these strands that have not separated from their junctions during the interval 
[t, t] . The amount 

dv 

- — (t,0,v) dt 

is the number of active strands (per unit mass) that detach (for the first time) from the network 
within the interval [t,t + dt], while the quantity 



d 2 u 

"dtd7 



(t, r, v) dtdr 



equals the number of strands (per unit mass) that have last merged with the network within 
the interval [t, r + dr] and separate from the network (for the first time after merging) during 
the interval [t, t + dt]. 

The rate of detachment <£> is defined as the ratio of the number of active strands that separate 
from temporary junctions per unit time to the total number of active strands. Applying this 
definition to active strands that were connected with the network at the initial instant t — 0, 
and to those that merged with the network within the interval [r, r + dr] , we arrive at the 
differential equations 



dv 

— (t,0,v) = -$(v)v(t,0,v), 



d 2 v , . ^, dv . , 
— (t,r,,) = -$(,)-(t,r,,). 



(31) 



Integration of Eq. (31) with initial conditions (29) (where we set t — 0) and (30) implies that 



v(t,0,v) = iYajo(i;)exp -$(v)t 



Ov r 

— (t,T,v) = <p(T,v)exp\-$(v)(t - t) 



To exclude the function <p(t,v) from Eq. (32), we use the identity 

/■* dv 

v(t, t, v) = u(t, 0, v) + J o —(t,T,v)dr. 
Substitution of expressions (29) and (32) into Eq. (33) results in 

N a p(v) — N a p(v) exp -$(v)t + / ip(r,v) exp -$(v)(t - t) 



dr. 



(32) 



(33) 



(34) 



The solution of linear integral equation (34) reads <p(t,v) = N a p(v)^(v). It follows from this 
equality and Eq. (32) that 



dv 



(t,r,v) = N a p(v)$(v)exp —$(v)(t - t) 



(35) 



We adopt the conventional assumptions that (i) the excluded-volume effect and other multi- 
chain effects are screened for individual strands by surrounding macromolecules, (ii) the energy 
of interaction between strands can be taken into account with the help of the incompress- 
ibility condition, and (iii) thermal oscillations of junctions can be disregarded, and the strain 
tensor for the motion of junctions at the micro- level coincides with the strain tensor for macro- 
deformation. 
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At isothermal deformation with small strains, a strand is treated as an isotropic incom- 
pressible medium. The strain energy of an active strand Wq is determined by the conventional 
formula 

Wq = pe : e, 

where p is an average elastic modulus of a strand, e is the strain tensor for transition from 
the reference (stress-free) state of the strand to its deformed state, the prime stands for the 
deviatoric component of a tensor, and the colon denotes convolution of two tensors. 

According to the affinity hypothesis, the strain energy w (t, 0) of an active strand that has 
not separated from the network during the interval [0, t] reads 

where e(t) is the strain tensor for transition from the initial (stress- free) state of the network to 
its deformed state at time t. With reference to [21] , we suppose that stress in a dangling strand 
totally relaxes before this strand captures a new junction. This implies that the stress-free state 
of an active strand that merges with the network at time r > coincides with the deformed 
state of the network at that instant. The mechanical energy of an active strand that has last 
merged with the network at time r e [0, t] is given by 

w{t,r) = /i[e(t)-e(r)]': [e(t) - e(r) 

Multiplying the strain energy per strand by the number of active strands per unit mass and 
summing the mechanical energies of active strands linked to temporary junctions with various 
activation energies, we find the strain energy per unit mass of an equivalent network 



W(t) = fi \u(t,0,v)e'(t):e'(t) + 



* dv 
o dr 



(t,r,v) e(t)-e(r) : e{t) - e(r) 



dr >dv 



(36) 



Differentiating Eq. (36) with respect to time t and using Eqs. (32), (33) and (35), we arrive at 
the formula 

,IW (t) = A(t):^-(t)-B(t), (37) 



dt 



dt 



where 



A(t) 

B{t) 



2fiNAe(t)- / e(r)dr / $(w)exp -$(«)(*- r) p(v)dv 



— (t,T,v)[e(t)-e(T)\ : [e(t) 



e t 



dr \dv > 0. 



(3* 



(39) 



For isothermal deformation of an incompressible medium, the Clausius-Duhem inequality reads 

Q 



dW V_ _ d£_ > 
dt p dt ~ 



where p is density, Q is internal dissipation per unit mass, and a stands for the stress tensor. 
Substitution of Eq. (37) into this equation implies that 



Q(t) = -[a\t)-pA{t) 



d£_ 
~dl 



(t)+B(t) > 0. 



(40) 
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As the function B(t) is non-negative, see Eq. (39), dissipation inequality (40) is satisfied, 
provided that the expression in the square brackets vanishes. This assertion together with Eq. 
(38) results in the constitutive equation 



a(t) = -P(t)I + 2G\ e(t) - / e\r)dr / exp - r) p(v)dv 



(41) 



where P(t) is pressure, I is the unit tensor, and G = p/iiV a is an analog of the shear modulus. 
Formula (41) describes the time-dependent response of an equivalent network at arbitrary 
three-dimensional deformations with small strains. In what follows, we confine ourselves to 
shear tests with 

e(t) = e(t)eie 2 , 

where e(t) is the shear strain, and e m (m = 1,2,3) are unit vectors of a Cartesian frame. 
According to Eq. (41), the shear stress a(t) is given by 

(42) 



a(t) = 2G|e(t) - J e(r)rfr J exp[-$(v)(t - T)\p(v)dv }. 

It follows from Eq. (42) that in a shear oscillation test with 

e(t) = e exp(iujt), 

where e and u are the amplitude and frequency of oscillations, and i 
complex modulus 



-1, the transient 



is determined by the formula 



2e(t) 



&(v)p(v)dv / exp 



$(v) + iuj)s 



ds 



where s = t — r. This equality implies that the steady-state complex modulus 



G*(u) = lim G*(t,u) 

t— >oo 



is given by 



G*(u) =G 



-p(v)dv. 



$(w) + IUJ" 

This equality together with Eq. (28) implies that the steady-state storage G'iui) and loss G"(uS) 
shear moduli read 

,2 



G"(u) 



G 
G 



uj- 



o $gexp(-2'u) + uj 2 

$o exp(— v)uj 
o $gexp(-2f) + cj 2 



p(v)dv, 
p(v)dv. 



(43) 



To fit the experimental data, we adopt the random energy model jUj with the quasi-Gaussian 
distribution function p(v), 



p(v) = po exp 



{v - Vf 



2S 2 



(v > 0), 
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p(v) = (v < 0), 



(44) 



where V and E are adjustable parameters (the apparent average activation energy and the 
apparent standard deviation of activation energies, respectively), and the constant p is found 
from the normalization condition 

p{y)dv = 1. (45) 



/ 

Jo 



Governing equations (43) and (44) involve four material constants: (i) the instantaneous shear 
modulus G, (ii) the attempt rate for rearrangement of strands $o, (hi) an analog of the average 
activation energy for rearrangement of strands in a network V, and (iv) an analog of the 
standard deviation of activation energies E. 

When the dimensionless ratio £ = E/V is small compared to unity (it will be shown later 
that this condition is satisfied for our experimental data), the number of adjustable parameters 
may be reduced to three. Assuming that 



V 



«1, 



(46) 



we can employ the first equality in Eq. (44) for an arbitrary (positive and negative) v. Replacing 
the lower limit of integration in Eqs. (43) by — oo, we obtain 



G» = G Po J 
G» = G Po 



-oo $0 exp(-2w) + u 2 
30 $ exp(—v)u 
oo $g exp(-2v) + oj 2 



exp 



cxp 



(v - V f 



2E 2 
(t; - V) 2 



2E 2 



dv, 
dv, 



(47) 



where 



1 



Po 



27TVS 

To exclude the attempt rate <3> from the consideration, we introduce the notation 



$ = $*exp(v ), 

where is a given value (in what follows, we set = 10 10 s _1 ), and v 
stituting this expression into Eq. (47) and introducing the new variable v' 
that 



ln$ /$^. Sub- 
v — v , we find 



G"{u) 



Gp 
Gp 



oo $ 2 exp(-2w') + to 2 
oo $lexp(-2v') + lo'- 



cxp 



■ cxp 



V 



f\2 



2E 2 
(v' - V') 2 



2S 2 



dv', 
dv', 



where 



V' = V - v = V 



ln^. 



Omitting primes for the sake of simplicity and replacing the lower limits of integration by zero, 
we return to Eqs. (43), where the unknown attempt rate $ is replaced by $*. This implies 
that each set of observations for the storage and loss shear moduli, G'(uo) and G"(uo), is entirely 
determined by three quantities: G, V and S. For hybrid nanocomposites subjected to thermal 
treatment, these parameters are functions of annealing temperature T a and annealing time t a . 
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Fitting of observations 



To assess the effect of temperature and time of annealing, we determine the quantities G, V and 
£ by matching the experimental data depicted in Figures 1, 2, 4, 5, 7, 8, 10 and 11. Each set 
of observations for G'(uo) and G"(u) is approximated separately. We fix some intervals [0, Knax] 
and [0, E max ], where the "best-fit" parameters V and £ are assumed to be located, and divide 
these intervals into J subintervals by the points = iAV and = jAE (i, j = 1, . . . , J— I) 
with Ay = V max /J and AS = £ max /J. For any pair {V {i \ E^}, the coefficient p in Eq. (44) 
is calculated from Eq. (45), where the integral is evaluated numerically by Simpson's method 
with 400 points and the step At> = 0.1. The integrals in Eq. (43) are calculated by using the 
same technique. The shear modulus G is found by the least-squares method from the condition 
of minimum of the function 



| [^cxp(^m) - G' num (u; m ) + Gg Xp (u; m ) - G" um (u; m ) j 



where the sum is calculated over all experimental points u> m , G' and G" xp are the storage 
and loss moduli measured in a test, and G' num and G^ um are given by Eq. (43). The "best-fit" 
parameters V and E are determined from the condition of minimum of the function F on the 
set {\/ (i) ,E^ = 1,...,J- 1)}. After finding the "best-fit" values V {i) and E^, this 

procedure is repeated twice for the new intervals [V^, V^] and [E^" 1 ), £^ +1 )], to ensure 
an acceptable accuracy of fitting. Figures 1, 2, 4, 5, 7, 8, 10 and 11 demonstrate excellent 
agreement between the experimental data and the results of numerical simulation. 

For each annealing time t a and annealing temperature T a , we (i) find the ratio of mass- 
average molecular weights d w from Eqs. (1) and (2) and the observations depicted in Figures 3, 
6, 9 and 12, (ii) calculate the material constants G, V and E by matching the experimental data 
reported in Figures 1, 2, 4, 5, 7, 8, 10 and 11, and (iii) plot the quantities G, V and E versus 
d w in Figures 17 and 18. The experimental data are approximated by the linear equations 

G = G + Gid w , V = V , E = E + S 1 rf w , (48) 

where the coefficients G i: Vi and Ej (i = 0, 1) are calculated by the least-squares method. 
Figures 17 and 18 show that Eq. (48) correctly describes changes in the adjustable parameters 
with mass-average molecular weight. The average activation energy for separation of strands 
from temporary junctions V is independent of molecular weight, whereas the shear modulus G 
and the standard deviation of activation energies E noticeably grow with M w . It can also be 
seen from Figure 18 that inequality (46) is satisfied with a reasonable level of accuracy (the 
ratio on left-hand side of Eq. (46) does not exceed 0.2). 



Discussion 

We begin with the analysis of material constants in Eq. (1) that describes the effect of frequency 
of oscillations u> on the modulus of complex viscosity rj. According to Figure 15, the exponent 
a is practically independent of mass-average molecular weight. It follows from this conclusion 
and Eq. (1) that the curves rj{uj) measured at various annealing temperatures T a and annealing 
times t a and plotted in double-logarithmic coordinates may be superposed (with a high level of 
accuracy) by shifts along the horizontal (frequency) and vertical (viscosity) axes. 
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Figure 16 demonstrates that the characteristic time r linearly decreases with mass-average 
molecular weight. This implies that the reciprocal quantity r _1 (with the dimension of rate) 
linearly grows with M w . The latter conclusion is in agreement with observations by Bywater 
and Black [45J, who found a similar trend for the degradation rate of poly(methyl methacrylate) 
and poly(a-methylstyrene). 

Figure 17 shows that the elastic modulus G linearly increases with mass-average molecular 
weight. A linear relation between the shear modulus and the molecular weight between entan- 
glements M e provides a basis for the statistical theory of rubber elasticity. It follows from this 
dependence and Figure 17 that the molecular weight between entanglements M e is proportional 
to the mass-average molecular weight M w . The latter result appears to be quite natural. 

According to Figure 18, the average energy for detachment of strands from temporary 
junctions V is independent of molecular weight. This result seems natural as well. Indeed, for 
a melt of a hybrid nanocomposite, the average activation energy V may be treated as the energy 
of thermal fluctuations necessary for mutual displacement of two chains (or for detachment of 
a chain from a stack of clay platelets) to a distance at which the two chains (or the chain and 
the filler particle) weakly affect each other. It seems plausible to assume that this parameter is 
independent of the chains' length (because it reflects local interactions between their segments). 
The latter implies that thermal degradation of a nanocomposite melt (modelled as scission of 
chains and annihilation of end- and side groups) should not affect V, which is confirmed by the 
experimental data presented in Figure 18. 

Figure 18 reveals that the standard deviation of activation energies E linearly increases with 
mass-average molecular weight. To provide an explanation for this observation, we recall that 
the standard deviation of activation energies E may be thought of as a measure of heterogeneity 
of an equivalent network of macromolecules, see Eq. (44). Thermal degradation of a polymer 
matrix results in homogenization of the network (as the rates of scission of macromolecules and 
detachment of side-groups are proportional to chains' lengths), which is observed as a reduction 
in E with a decrease in d w . It is worth noting that a similar decrease in the inhomogeneity of 
an equivalent network driven by thermal degradation of neat iPP was previously observed as a 
reduction in the polydispersity index with annealing time 

Figures 15 to 18 show that the adjustable parameters in the model are not affected by the 
history of thermal pre-treatment, but are determined by the current mass-average molecular 
weight M w exclusively. Some scatter should, however, be mentioned of the experimental data 
depicted in these figures. It may be explained by the fact that at each temperature T a and each 
annealing time t a , a new sample was used for testing, whose physical properties do not exactly 
coincide with those of other specimens. 

Concluding remarks 

A series of torsional oscillation tests with small strains have been performed at the temperature 
T = 230 °C on a hybrid nanocomposite with a polypropylene matrix reinforced with 5 wt.% of 
MMT clay. Prior to mechanical tests, specimens were annealed at the temperatures T a = 250, 
270, 290 and 310 °C for various amounts of time t a ranging from 15 to 420 min. Thermal 
treatment induced thermal degradation of samples observed as a pronounced decrease in their 
mass-average molecular weight with exposure time t a . 
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With reference to the fragmentation-annihilation concept, a kinetic equation has been devel- 
oped for the concentration of chains with various lengths. This relation involves two adjustable 
parameters that are found by matching the experimental data for the evolution of mass-average 
molecular weight. In addition to the numerical analysis, an explicit solution of the kinetic equa- 
tion has been derived. 

A constitutive model has been developed for the viscoelastic response of a nanocomposite 
melt at isothermal three-dimensional deformations with small strains. The melt is treated as 
an equivalent transient network of strands bridged by temporary junctions. Its time-dependent 
behavior is modelled as separation of active strands from their junctions and attachment of 
dangling strands to the network. The rearrangement events occur at random times, when 
appropriate strands are thermally activated. 

Stress-strain relations for an equivalent heterogeneous network of strands (where different 
junctions have different activation energies for rearrangement of strands) have been derived by 
using the laws of thermodynamics. The constitutive equations involve three material parame- 
ters that are determined by matching the experimental data for the storage and loss moduli as 
functions of frequency of oscillations uj. Fair agreement is demonstrated between the observa- 
tions and the results of numerical simulation. 

The following conclusions are drawn: 

1. The average activation energy for rearrangement of strands in a transient network V is 
practically independent of molecular weight. 

2. The standard deviation of activation energies S and the shear modulus G linearly grow 
with mass- average molecular weight M w . 

3. The values of these parameters are independent of the history of thermal pre-treatment 
and are entirely determined by the current mass-average molecular weights. This implies 
that thermal degradation may be used as a quick-and-dirty test directed to establish 
correlations between the viscoelastic properties and molecular weights of polymers and 
hybrid nanocomposites. 
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List of figures 

Figure 1: The storage modulus G' versus frequency u. Circles: experimental data on iPP 
MMT nanocomposite annealed at T a = 250 °C for t a = 0, 60 and 90 min, from top to 
bottom, respectively. Solid lines: results on numerical simulation 

Figure 2: The loss modulus G" versus frequency u. Circles: experimental data on iPP-MMT 
nanocomposite annealed at T a = 250 °C for t a = 0, 60 and 90 min, from top to bottom, 
respectively. Solid lines: results on numerical simulation 

Figure 3: The complex viscosity rj versus frequency uo. Circles: experimental data on iPP- 
MMT nanocomposite annealed at T a = 250 °C for t a = 0, 60 and 90 min, from top to 
bottom, respectively. Solid lines: results on numerical simulation 

Figure 4: The storage modulus G' versus frequency oo. Circles: experimental data on iPP 
MMT nanocomposite annealed at T a = 270 °C for t a = 0, 60 and 90 min. Solid lines: 
results on numerical simulation 

Figure 5: The loss modulus G" versus frequency uo. Circles: experimental data on iPP-MMT 
nanocomposite annealed at T a = 270 °C for t a = 0, 60 and 90 min. Solid lines: results on 
numerical simulation 

Figure 6: The complex viscosity rj versus frequency uo. Circles: experimental data on iPP 
MMT annealed at T a = 270 °C for t a = 0, 60 and 90 min, from top to bottom, respectively. 
Solid lines: results on numerical simulation 

Figure 7: The storage modulus G' versus frequency uo. Circles: experimental data on iPP 
MMT nanocomposite annealed at T a = 290 °C for t a = 0, 15, 30, 60, 120 and 420 min, 
from top to bottom, respectively. Solid lines: results on numerical simulation 

Figure 8: The loss modulus G" versus frequency uo. Circles: experimental data on iPP-MMT 
nanocomposite annealed at T a = 290 °C for t a = 0, 15, 30, 60, 120 and 420 min, from top 
to bottom, respectively Solid lines: results on numerical simulation 

Figure 9: The complex viscosity 77 versus frequency uo. Circles: experimental data on iPP 
MMT nanocomposite annealed at T a = 290 °C for t a = 0, 15, 30, 60, 120 and 420 min, 
from top to bottom, respectively. Solid lines: results on numerical simulation 

Figure 10: The storage modulus G' versus frequency uo. Circles: experimental data on iPP 
MMT nanocomposite annealed at T a = 310 °C for t a = 0, 15, 30 and 60 min, from top to 
bottom, respectively. Solid lines: results on numerical simulation 

Figure 11: The loss modulus G" versus frequency uo. Circles: experimental data on iPP 
MMT nanocomposite annealed at T a = 310 °C for t a = 0, 15, 30 and 60 min, from top to 
bottom, respectively. Solid lines: results on numerical simulation 
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Figure 12: The complex viscosity 77 versus frequency u. Circles: experimental data on iPP 
MMT nanocomposite annealed at T a = 310 °C for t a = 0, 15, 30 and 60 min, from top to 
bottom, respectively. Solid lines: results on numerical simulation 

Figure 13: The ratio of mass-average molecular weights d w versus annealing time t a . Symbols: 
experimental data on iPP-MMT nanocomposite annealed at temperatures T a °C. Unfilled 
circles: T a = 250; filled circles: T a = 270; asterisks: T a = 290; stars: T a = 310. Solid line: 
results of numerical simulation 

Figure 14: The shift factor A versus annealing temperature T a . Circles: experimental data 
on iPP-MMT nanocomposite. Solid line: approximation of the experimental data by Eq. 
(3) with A = 13.32 and A x = 7 A3 ■ 10 3 

Figure 15: The dimensionless exponent a versus the ratio of mass-average molecular weights 
d w . Symbols: treatment of observations on iPP-MMT nanocomposite at annealing tem- 
peratures T a °C. Unfilled circles: T a = 250; filled circles: T a = 270; asterisks: T a = 290; 
stars: T a = 310. Solid line: approximation of the experimental data by Eq. (4) with 
«o = 0.68 and a\ = 0.17 

Figure 16: The characteristic time r versus the ratio of mass-average molecular weights d w . 
Symbols: treatment of observations on iPP-MMT nanocomposite at annealing tempera- 
tures T a °C. Unfilled circles: T a = 250; filled circles: T a = 270; asterisks: T a = 290; stars: 
T a = 310. Solid line: approximation of the experimental data by Eq. (4) with r = —0.39 
and T\ = 0.95 

Figure 17: The instantaneous shear modulus G versus the ratio of mass-average molecular 
weights d w . Symbols: treatment of observations on iPP-MMT nanocomposite at anneal- 
ing temperatures T a °C. Unfilled circles: T a = 250; filled circles: T a = 270; asterisks: 
T a = 290; stars: T a = 310. Solid line: approximation of the experimental data by Eq. 
(48) with G = -0.07 and d = 0.32 

Figure 18: The average activation energy for rearrangement of strands V and the standard 
deviation of activation energies £ versus the ratio of mass-average molecular weights d w . 
Symbols: treatment of observations on iPP-MMT nanocomposite at annealing temper- 
atures T a °C. Unfilled circles: T a = 250; filled circles: T a = 270; asterisks: T a = 290; 
stars: T a = 310. Solid lines: approximation of the experimental data by Eq. (48) with 
V = 15.82 (curve 1) and S = 0.95 and Si = 1.93 (curve 2) 
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